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Abstract 

The computation of the 0{as) QCD corrections to the cross sections for W^jj 
production in hadronic collisions is presented. We consider the case of real photons 
in the final state, but include full leptonic decays of the W. Numerical results for the 
LHC and the Tevatron are obtained through a parton level Monte Carlo based on the 
structure of the VBFNLO program, allowing an easy implementation of general cuts 
and distributions. We show the dependence on scale variations of the integrated cross 
sections and provide evidence of the fact that NLO QCD corrections strongly modify 
the LO predictions for observables at the LHC both in magnitude and in shape. 
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1 Introduction 



Precise and reliable predictions of cross sections at hadron colliders require the calcu- 
lation of higher order QCD corrections. As part of such a program we have in the past 
determined next-to-leading-order (NLO) QCD corrections for the production cross sections 
of various combinations of three electroweak bosons, including W^W~ Z [Ij, W^IV^W^ 
and W^ZZ [2], W^W'j^ZZj) and W^Z^ production [4J. In all cases, leptonic decays 
of the weak bosons were included in the calculations. For the production of three weak 
bosons, these results were verified against independent calculations [HE] which are available 
for on-shell weak bosons and neglecting Higgs boson exchange. 

In the present paper, we present the 0[a^ QCD corrections for processes with two 
photons and one W in the final state, namely the production of 

"iy*"77" pp, pp ^ uil'^'yy +X 

"W^"77" pp, pp ^ l^ui'yy +X. (1.1) 

A calculation of Vr+77 production at NLO QCD accuracy has been performed before [7], 
including the <? — ?■ 5'7 fragmentation contributions. However, also Ref. [7] treats the W as 
stable. By including off shell effects and, in particular, photon radiation from the final state 
charged lepton we aim at providing complementary information. 

Processes involving two or three electroweak bosons in the final state are relevant for 
studying anomalous gauge interactions, as they give direct access to triple and quartic cou- 
plings [8j. Vr^77 production is sensitive to the WWj and WW'j'j vertices. In addition, a 
final state with two photons and missing transverse energy is relevant in a variety of beyond 
the standard model scenarios |H]: in gauge- mediated supersymmetry breaking, for instance, 
the neutralino is often the next-to-lightest supersymmetric particle and decays into a photon 
plus a gravitino, giving a signal of two photons and missing Et- In Ref. [10], a study of 
the background estimates for supersymmetry motivated di-photon production searches has 
been performed, pointing out the relevance of the W^'yy production process as a standard 
model (SM) background in case of electron misidentification. Another possible application 
is an estimate of backgrounds when searching for WH production, followed by Higgs decay 
to two photons. 

The present work closely follows our previous calculations of NLO QCD corrections to 
triple electroweak boson production. In particular, photon isolation is implemented within 
the Frixione approach [11] and we, thus, avoid the need for the inclusion of g — )■ g7 fragmen- 
tation contributions, which were discussed for IV "'"77 production by Baur et al. Similar 
to previous work on triple weak boson production, we find that the QCD corrections are 
sizable and also modify the shape of the differential distributions for many observables: this 
proves that a simple rescaling of the LO results is not adequate and a full NLO Monte Carlo 
is needed for any quantitative determination of quartic couplings at the LHC. We have im- 
plemented our calculation within the VBFNLO framework {121, a parton level Monte Carlo 
program which allows the definition of general acceptance cuts and distributions. 

The paper is organized as follows: in Section [2] we show the Feynman diagrams relevant 
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for the calculation and provide an overview of the strategies used to compute the real and 
virtual corrections and the various checks performed to ensure the numerical accuracy of 
our code. In Section [31 we show numerical results, including the scale variations of the LO 
and NLO integrated cross sections and some selected differential distributions. Particular 
concern will be given to the consequences of approximate radiation zeroes which largely 
disappear when going from LO to NLO distributions, as was already emphasized in Ref. [7]. 
Conclusions are given in Section HI 




For the calculation of virtual corrections, it is convenient to classify the contributing 
Feynman graphs according to the number of electroweak boson vertexes which are attached 
to the quark lines: graphs with three vertexes lead to pentagon loops, while boxes and 
triangles are the most complex structures for graphs with two or one electroweak vertex on 
the quark line. Examples for the three classes of diagrams contributing to W^'yy production 
at Born level are shown in Fig. [TJ Diagrams (a) and (b) exhibit the dependence on the 
quartic (WW^yj) and triple (14^1^7) gauge boson couplings. 

As is customary in all VBFNLO calculations, we include the full spin correlations of the 
leptons coming from the W decay and also the final state photon radiation off the charged 
leptons. As in all previous calculations, we have factorized the subgraphs involving boson 
splitting and decay into leptonic tensors [13], computed once per phase-space point in order 
to speed up the code, and we made use of the helicity technique introduced in Ref. [11] for 
the computation of matrix elements. The cancellation of infrared divergences coming from 
real and virtual corrections at NLO is achieved using the Catani-Seymour dipole subtraction 
method [15]. We note here that the strategy of separately computing the leptonic tensors 
proves particularly useful given the large number (372) of real emission diagrams at NLO. 

For the real emission contributions, there is an additional issue to consider: the infrared 
singularity that may result from a photon emitted collinearly to a massless quark. The simple 
rejection of any event containing partons within a cone drawn around the photon direction 
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would inevitably spoil the cancellation of infrared divergences for subprocesses with a gluon 
in the final state. As a solution, we use the photon isolation cut proposed by Frixione jlT] 
which is described in more details in Section [3J 

The NLO virtual corrections result from one-loop diagrams obtained by attaching a gluon 
line to the quark-antiquark line in diagrams like the ones depicted in Fig. [H We combine the 
virtual corrections into three different groups, which include all loop diagrams derived from 
a given Born level configuration [13j. This leaves us with three universal building blocks, 
namely corrections to one, two or three vector bosons attached to the quark line. In all 
cases, the infrared divergent contributions factorize in terms of their corresponding Born 
amplitude and the full virtual term is given by 



Mv = Mv + 



ra + e) 



2 



3 
e 



+ 



4 it' 



(2.1) 



where Q is the partonic center-of-mass energy, i.e^the invariant mass of the W^jj system, 
Mb the total born contribution, and the term My includes the finite parts of the virtual 
corrections to 2 and 3 weak boson amplitudes which we call virtual-box and virtual-pentagons 
in the following, so named since boxes and pentagons constitute the most complex loop 
diagrams, respectively. These finite terms can be calculated in d = 4 dimensions. The tensor 
coefficients of the loop integrals have been computed by means of the Passarino-Veltman 
reduction formalism up to the box level, but avoiding the explicit calculation of inverse 
Gram matrices by solving a system of linear equations which is a more stable procedure close 
to singular points. For pentagons, we use the Denner-Dittmaier reduction formalism [T6] . 
which is now fully implemented in the public version of the VBFNLO code for all multi- 
boson processes. As in our previous calculation |3], we explicitly checked that there are no 
additional infrared singularities other than those proportional to the Born amplitude. 

A powerful test of the virtual corrections is provided by Ward identities which connect 
virtual-pentagon and virtual-box contributions. For this purpose we write the polarization 
vector of the W as 



(2.2) 



where is the momentum of the W and the remainder e^y is chosen in such a way that 



€w ■ (^71 +^72) = 



(2.3) 



i.e., the time component of the shifted polarization vector is zero in the center-of-mass system 
of the photon pair. Via Ward identities, a pentagon contracted with the W momentum, q^, 
can always be reduced to a difference of boxes: in this way we "shift" a fraction of pentagon 
diagrams to box subroutines. The fact that the sum of the virtual contributions does not 
change upon this shift provides a very powerful consistency check of our implementation. 

The numerical accuracy of our code for tree level amplitudes has been tested against 
MadGraph [TT] at the level of amplitudes and against Sherpa [18] for integrated cross sections, 
finding agreement at the level of machine precision and at the per mill level, respectively. 
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As a final test, we have made a comparison with the numbers of the proceeding paper of 
Ref. [7] for the production of on-shell bosons without leptonic decays. For that we have 
neglected all diagrams with the photon emitted from the lepton line and used narrow-width 
approximation for the vector boson decay. In Table [H we show the comparison for W'^'-ff+X 
production at LO at the LHC We do not present an explicit comparison at NLO since we 
do not include fragmentation functions and our isolation procedure for the photon differs 
from the one used in Ref. [?|. As input parameters for this comparison, we use a at = 
and standard FormCalc EW parameters. 



Scale 


Program 


LO [fb] 


100 GeV 


VBFNLO 
Ref. trj 


7.317(1) 
7.253(5) 



Table 1: Comparison between our results and the ones of Ref . [7j for pp — )■ W'^'j'j +X at 
LO at the LHC. The input parameters and settings are taken from Ref. [7]. 



Additionally, to control the numerical stabihty of our code, a gauge test based on Ward 
identities has been used throughout the virtual implementation for each phase space point. 

3 Results 

3.1 Definition of scales and cuts 

We have implemented our calculation as a NLO Monte Carlo program based on the structure 
of the VBFNLO code [12] . For the electroweak parameters, we use the W and Z boson masses 
and the Fermi constant as input. From these, we derive the electromagnetic coupling and 
the weak mixing angle via tree level relations, i.e. we use 

mw = 80.398 GeV mz = 91.1876 GeV sin^ (Ow) = 0.22264 

Gf = 1.16637 ■ 10"^ GeV-2 a'^ = 132.3407. (3.1) 

We do not consider bottom and top quark effects. The remaining quarks are assumed 
to be massless and we work in the approximation where the CKM matrix is the identity 
matrix. We choose the invariant "iy^77" mass as the central value for the factorization and 
renormalization scales: 

/iF = /^i? = /io = \/ {Pi + Pu+P-ri +^72)^- (3-2) 

We use the CTEQ6L1 [19] parton distribution function at LO and the CTIO [20] set with 
as{mz) = 0.1180 at NLO. 

A real photon in the final state can be emitted either from the initial quark line or 
from the final-state charged lepton. For efficient Monte Carlo generation, we divide the 
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phase space into three separate regions to consider all the possibilities and then sum the 
contributions to get the total result. The regions are generated as triple electroweak boson 
production as well as Wji and W'j2 production with (approximately) on-shell — )■ I'^i'ij 
(or W~ — )■ three-body decay, respectively. 

We impose a set of minimal cuts on the rapidity, , and the transverse momenta, Pt£(7) , 
of the charged leptons and the photon, which are designed to represent typical experimental 
requirements. Furthermore, leptons, photons and jets must be well separated in the rapidity- 
azimuthal angle plane. Specifically, the cuts imposed are 

Pt^(^) > 20 GeV |y£(^)|<2.5 i?^^ > 0.4 i?^^ > 0.4 > 0.4 Rj^ > 0.7 

(3.3) 

where, in our simulations, a jet is defined as a colored parton of transverse momentum 
PTj > 30 GeV and rapidity \yj\ < 4.5. 

Since we do not include any fragmentation contribution, we must find a way to reject 
events in which a quark and a photon are collinear in the final state, i.e. we have to provide 
a prescription to isolate the photons. We choose to implement the procedure defined in 
Ref. [11] : if i is a parton with transverse energy Et^ and has a separation Ri^ with a photon 
of transverse momentum px-y, then the event is accepted only if 

1 — cos 

Et, 9{6 - Ri^) < (for all 5 < 5o) (3.4) 

1 — cos do 

where 6o is a fixed separation which we set equal to 0.7. A quick look at Eq. (13. 4 p reveals that 
a sufficiently soft parton can be arbitrarily close to the photon axis, while the energy of an 
exactly collinear parton must be vanishing in order to pass the isolation cut. Collinear-only 
events (leading to fragmentation contributions) are thus rejected while soft emissions are 
retained as desired. 



3.2 Integrated results 

In Table [21 we give results for the integrated cross sections for W'^'-fy production at the LHC 
(14 TeV) for the given cuts as well as for a harder cut on the photon transverse momentum 
of pT'y > 30 GeV: the NLO corrections are large, enhancing the LO result by more than a 
factor of 3 in all cases as can be seen by the K-factor defined as K= a^^^ /a^'^ . This value 
is larger than the ones in other triple vector boson production channels, where we observe 
K-factors between 1.5 and 2. However, it is consistent with the results of Ref. [7], where a 
K-factor of 2.93 was found after isolation cuts (which differ from ours). We already note 
here that the large K-factor can be explained by the suppression of the LO cross section due 
to the so called radiation zero. This will be further investigated in Section 

In Table |3l the numbers for the Tevatron (1.96 TeV) are presented for the cuts as in 
Eq. (13. 3p (and for photon isolation as in Eq. (I3.4p ). but for less restrictive transverse mo- 
mentum cuts, Pt'^^i) > 10(10) GeV and pr-yii) > 20(10) GeV. We only show the W^'j'j case, 
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T HP L /e — 14 TpA/^I 
ijnv^ IV* — lev J 


LiVJ [iUJ 


NT n [fhi 


lY-icLClUi 


(t('W+77" e+z/e77) 








> 20(20) GeV 


2.529 


7.940 


3.14 


> 30(20) GeV 


0.979 


3.172 


3.24 


cr("iy"77" e~z/e77) 








> 20(20) GeV 


1.946 


6.759 


3.47 


> 30(20) GeV 


0.686 


2.583 


3.77 



Table 2: Total cross sections at the LHC for pp — )■ l^=''77+X with leptonic decays, at LO 
and NLO, and for two sets of cuts. Relative statistical errors of the Monte Carlo are below 

lo-l 



Tevatron 


LO [fb] 


NLO [fb] 


K-factor 


a('W+77" ^ e+z/e77) 








PT-yii) > 10(10) GeV 


4.779 


7.558 


1.58 


Pr7W > 20(10) GeV 


0.5591 


0.9415 


1.68 



Table 3: Total cross sections at the Tevatron for pp — W^'y^y + X with leptonic decays, at 
LO and NLO, and for two sets of cuts. Relative statistical errors of the Monte Carlo are 
below 10~^. 



because W 77 production is exactly symmetrical at a pp collider. The NLO enhancement 
is less pronounced for this collider, but still amounts to 60-70% of the LO result. 

In the following, we include a combinatorial factor of 2 in all figures which corresponds 
to the production of electrons and muons. We have studied the scale uncertainty of the total 
cross section by varying the renormalization and factorization scales as 

fiF,m = ^-f^o (0.1<e<10). (3.5) 

In Fig. [2]and Fig. [3], we show numerical results for W'^'-ff production and W^jj production, 
respectively, within the cuts of Eqs. (13. 3[ 13. 4p . On the left panel of each figure, we show 
the overall scale variation of our numerical predictions at LO and NLO: as already shown 
in Table [21 the NLO K-factor is large both in absolute value (~ 3) and compared to the 
LO scale variation. The NLO scale uncertainty is about 10% when varying the factorization 
and the renormalization scale /i = /ii? = /i^ up and down by a factor 2 around the reference 
scale fj,o = mwyy and is mainly driven by the dependence on fin, which gives a negative 
slope with increasing energy, while the dependence on shows the opposite behaviour. In 
the left panels, we also show results for additional jet veto cuts, requiring p^j < 50 GeV 
or pTj < 30 GeV. From these curves, we can see that a large contribution to the total 
cross section is due to real jet radiation. While it is evident that the renormalization scale 
variation is highly reduced by a jet veto, this reduction should not be interpreted as a smaller 
uncertainty of the vetoed cross section: a similar effect in WZj production could be traced 
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Figure 2: Left: Scale dependence of the total LHC cross section for pp — )■ iy~77 + X — )■ 
^~77 + ;^rp + X at LO and NLO within the cuts of Eqs. ( 13. 31 \3.4\) . The factorization and 
renormahzation scales are together or independently varied in the range from 0.1-/io to IO-hq. 
Right: Same as in the left panel but for the different NLO contributions at fJ.F = fJ-R = ^/^o- 



to cancellations between different regions of phase space and, thus, the small variation is 
cut-dependent [2T] . 
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Figure 3: Same as Fig. El but for pp 1^+77 + X £+77 + + X at the LHC. 

Since the 1-jet contributions to the 0{as) cross section are determined at LO only, 
their scale variation is large. In fact, most of the scale variation of the total NLO result 
is accounted for by the real emission contributions, defined here as the real emission cross 
section minus the Catani-Seymour subtraction terms plus the finite collinear terms. This 
is more visible in the right panels, where we show the scale dependence and compare the 
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Figure 4: Left: Invariant mass distribution of the photon pair for pp — )■ W'^'-ff + X — )■ 
i+^yy -\- -\- X production at the LHC. LO and NLO results are shown for /xf = yUj? = /^o 
and the cuts of Eqs. ^3.3\ \3.4\) . Right: K-factor as defined in Eq. / I3.6)) . 

size of the different parts of the NLO calculation. As for the relative size of the NLO 
terms, the real emission contributions dominate and are even larger than the LO terms plus 
virtual terms proportional to the Born amplitude. Non-trivial virtual contributions, namely 
the interference of the Born amplitude with virtual-box and virtual-pentagon contributions, 
represent less than 1% of the total result and their scale dependence is basically flat. 



3.3 Differential cross sections 



Our numerical results show that the NLO corrections have a strong dependence on the phase 
space region under investigation. As a consequence, a simple rescaling of the LO results 
with a constant K-factor is not allowed. As practical examples, we plot several differential 
distributions at LO and NLO together with the associated K-factor, defined as 

where x denotes the considered observable. We usually show only one of the W'^^yj and 
W~'-f'-f distributions in the following since they share similar behaviors. Also, we include in 
all figures the results with an additional jet veto cut, requiring ptj < 50 GeV. 

In Fig. m we show the differential cross section as a function of the invariant mass of the 
photon pair, M^^ for W'^'yy production. The corresponding K-factors lie between 2.5 and 
3.5 in most of the phase space region. 

In Fig. 15]®, we show the transverse momentum distribution of the harder (softer) photon 
PT'y,max{PT'y,min) in W'^'-yj productiou. The K-factor is almost constant in the case of the 
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harder photon px distribution. For the px distribution of the softer photon, we get the 
same large K-factor at low values of transverse momenta, but with a decrease of the NLO 
enhancement above pt ~ 50 GeV. 




Figure 5: Left: Transverse-momentum distribution of the hardest photon in pp — )■ + 
X — )■ i'^'j'j+^j^+X production at the LHC. LO and NLO results are shown for = /x_r = /xo 
and the cuts of Eqs. i3.3[ \3.4]] . Right: K-factor as defined in Eq. ( 13.6)) . 




Figure 6: Left: Transverse-momentum distribution of the softest photon in pp W'^^y^y + 
X — > production at the LHC. LO and NLO results are shown for ^ip = yUi?, = /^o 

and the cuts of Eqs. ( 13.31 \3.4\) . Right: K-factor as defined in Eq. ( 13. 

Similar features emerge for distributions of the photon-lepton separations. In Fig. [71 
the minimum of i?;^^ and -Rj^j is shown for W~'~f'~f production: the K-factor shows a large 
phase-space dependence and reaches values of 3 to 4 when the photons are emitted close to 
the lepton. 
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Figure 7: Left: Distribution of the minimum between Ri^^ and Ri-^^ in pp — )■ iy~77 + X — )■ 
^~77 + "^rp + X production at the LHC. LO and NLO results are shown for /xf = /^_r = /^o 
and the cuts of Eqs. ^3.3\ \3.4\) . Right: K-factor as defined in Eq. / I3.6)) . 



Much larger values of the K-factor are possible in phase space regions which are difficult 
to access due to the LO kinematics. An example with a K-factor larger than 20 is shown 
in Fig. [HI where the transverse-momentum distribution of the lepton-photon-photon system 
is plotted. This large K-factor, indeed, is a feature shared by other triple vector boson 
production channels: the increase of the transverse-momentum of the observable electroweak 
system is compensated by a hard jet which appears first at NLO. 




Figure 8: Left: Transverse-momentum distribution of the £"77 system inpp — )■ W'^y^y+X — )• 
^~77 + "^rp + X production at the LHC. LO and NLO results are shown for np = = hq 
and the cuts of Eqs. ( I3.3[ \3.4\) . Right: K-factor as defined in Eq. ( 13. 6j) . 
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Figure 9: Left: Dependence of the NLO cross section on the couple of variables 
(min(Mi^j, M^^J, -R^i^a)- Same scales and cuts as in the previous diagrams. Right: K-factor 
as defined in Eq. 113.6]] . 




Figure 10: Left: Dependence of the NLO cross section on the couple of variables 
{mm{Ri^j^, Ri^^) , R^j^^^) . Same scales and cuts as in the previous diagrams. Right: K-factor 
as defined in Eq. ( 13. 6j) . 
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Large K-factors also appear in the two-dimensional distributions of Figs. [9]and[T0l where 
the dependence of the cross section on the pairs of variables (min(M;^-^, M^^j), .R^^^j) and 
(min(i?;^^, i?;^2)5 -^7172)5 respectively, are shown, together with the corresponding K-factors. 
It is evident that very large NLO corrections arise particularly in regions where the lepton and 
photons are emitted close to each other and for low invariant masses of the lepton-photons 
system. However, since these regions contribute very little to the total cross section, this 
does not affect the overall K-factor, which is similar to the one observed in the bulk of 
the corrections in the differential distributions of Figs. H] to [HI Note, however, that for the 
study of anomalous couplings, frequently extreme kinematics are selected. Therefore, large 
K-factors might appear. 

3.4 Radiation Zero 

As is known from the general theorem of Ref. ^22j, the SM amplitude for the process qQ — >■ 
H^='=77 vanishes for cos 6*^ = ±1/3 when the two photons are coUinear. Here, 6^ denotes 
the angle between the incoming quark and the W boson in the parton-center of mass frame. 
This radiation zero is not present in the gluon-induced channels, which enter in this process 
at NLO and are important at the LHC due to the steep rise of the gluon pdfs with smaller 
Feynman-x. The radiation zero only remains present when additional neutral (e.g. gluonic) 
radiation is coUinear to the photon. Hence, additional QCD emission, as part of the NLO 
contribution to iy^77 production, is expected to spoil the radiation zero, similar to the 
production process [23] . The resulting strong increase of the cross section near cos = ± 1 /3 
can explain the large total K-factor for the W^jj process. 

The radiation zero can be investigated, following Ref. [23], via rapidity difference dis- 
tributions of the W and the photon pair system, y^^ — yw, where yw is obtained from the 
W momenta reconstructed out of the lepton and the missing transverse momentum. In our 
discussion we follow Ref. |7j and present results only for W'^'yy + X production since similar 
features are observed for W~'-fy + X production. The new element in the present discussion 
is that final state photon radiation (off the charged lepton) is included in the calculation. 

In Fig. [m we plot the distribution in rapidity separation between the W and the photon 
pair, for photons lying in the same (cos^ > 0) or opposite (cos^ < 0) hemispheres in the 
laboratory frame. Our results at LO differ from the left-hand plot shown in Figure 3 of 
Ref. [7], where the W is considered as stable: our zero-rapidity dip for cos6' < is stronger 
than that for cos 9 > 0, while the opposite behaviour is observed in Ref. In both cases, the 
NLO corrections almost fill the dips, making an observation of a radiation zero at the LHC 
very difficult. In the K-factor plot of Fig. [TlTc). the origin of the large total K-factor (~3) 
for W^77 production is clearly visible. This can be compared to other triple vector boson 
production processes with typical K-factor values between 1.5 and 2 pT-'B]. The suppression 
in the central region at LO due to the radiation zero is spoiled by the extra jet emission at 
NLO, giving rise to large K-factors in the bulk of the corrections and therefore large total 
K-factors. This is opposite to Figs. [8lfT0| where large K-factors appear in marginal phase 
space regions. 
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yy/'Vw* yyy'Vw* Vy/'yvv* 



Figure 11: Distribution in rapidity separation between the W and the photon pair, with 
the photons in the same (cos 6 > 0) or opposite (cos 9 < 0) hemispheres, at LO (left), NLO 
(center) and K- factor (right). 




yi7"yw+ Vf/'yvi* yi7"yw+ 



Figure 12: Same plots as in Fig. lTl\with the additional cut on Mt,iu > 70 GeV. 



The difference observed at LO in Ref. [7] as compared to our results is due to the extra 
radiation of photons from the lepton hne. In Fig. [121 we try to suppress this contribution by 
restricting the phase space for photon radiation in W decay. Imposing a cut, Mt^iu > 70 GeV, 
on the transverse mass of the charged lepton and neutrino, which is close to the kinematical 
limit for on-shell Ws, strongly disfavors final state photon emission. This results in y^^ — t/w 
distributions with a more pronounced dip for cos 6' > 0, in agreement with Ref. [7j. Note 
that even larger K-factors now appear in the central region, reaching values close to 20 in 
regions of large cross section. 

For a better comparison with Reference [7j, we also show our numerical results with an 
additional jet veto cut, ptj < 50 GeV, in Fig. [131 with or without the cut on the invariant 
mass of the lepton-neutrino pair. Note also that, following Ref. |24], we use true rapidity 
and not pseudorapidity distributions, as was done in Ref. [7]. We have checked that our 
pseudorapidity plots, once we unphysically remove the extra radiation of the photon from 
the lepton line, are similar to those of Ref. [7]. It is interesting to see in the last plots. 
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Figure 13: Left plot: NLO plot of Fig. [IT] with the additional jet veto cut ptj < 50 GeV. 
Right plot: NLO plot of Fig. [H] with the additional cuts prj < 50 GeV and Mtau > 70 GeV. 

similar to Fig. 3 of Ref. [7], that once extra radiation is vetoed, a small dip is visible at 
NLO. This feature is particularly apparent in the right figure of Fig. [131 where final state 
photon radiation is also reduced. However, whether this dip should be interpreted as a sign 
of the radiation zero is questionable since it is more pronounced for photons in opposite 
hemispheres. As was pointed out in Ref. [25], signing the quark direction according to 
overall boost of the event along the beam axis efficiently lifts the proton-proton initial state 
degeneracy at the LHC Thereby, it reduces the distortion to the radiation zero by QCD 
activity and might aid in a more significant observation of the radiation zero. We have not 
investigated these questions in the present work. 



4 Conclusions 



We have calculated the NLO QCD corrections to the processes pp, pp — )■ W 'yy+X with 
full leptonic and radiative decays of the W boson. This process can be relevant as a back- 
ground for New Physics searches and as a signal for the measurement of triple and quartic 
couplings at the LHC. It also represents an irreducible background to WH production, when 
the Higgs boson decays into two photons. 

Our numerical results show that this process gets very large corrections at NLO, even 
when compared to other triple boson production processes. The exceptionally large total K- 
factor is due to a cancellation between different diagrams at LO, leading to an approximate 
radiation zero, an effect which is spoiled by the extra jet emission at NLO. In certain regions 
of phase space the cross section increases by a factor of 20 with respect to the LO calculation. 
Apart from the vicinity of the radiation zero, this strong enhancement happens when the 
recoil against a hard jet at NLO is kinematically favorable. This latter effect was already 
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observed in the NLO corrections to other triple boson production processes. However, these 
high transverse momentum effects do not significantly contribute to the bulk of the cross 
section corrections and do not affect the total K-factor. Nevertheless, for the study of 
anomalous couplings, extreme kinematics are typically selected. Therefore, large K-factors 
might appear and the NLO corrections have to be taken into account. 

IV^77 production at the LHC provides an additional example of a cross section whose 
theoretical errors at LO are substantially underestimated by considering scale variations 
only: the LO factorization scale variation is much smaller than the size of NLO corrections. 
Remaining NLO scale variations are at the 10% level for the integrated W^jj production 
cross section at the LHC when varying fi^ = /i^ = /iq by a factor of 2 around the reference 
scale /io = mww'y 

Given the size of the higher-order corrections and, in particular, their strong dependence 
on the observable and on different phase space regions under investigation, a fully-flexible 
NLO parton Monte Carlo for W^jj production is required to match the expected precision 
of the LHC measurements. We plan to incorporate this and other processes with a final 
state photon into the VBFNLO package in the near future. 
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